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ABSTRACT 

We compute the angular momentum, the spin parameter and the related distribution function for Dark Matter 
halos hosting a spiral galaxy. We base on scaling laws, inferred from observations, that link the properties 
of the galaxy to those of the host halo; we further assume that the Dark Matter has the same total specific 
angular momentum of the baryons. Our main results are: (i) we find that the gas component of the disk 
significantly contributes to the total angular momentum of the system; (ii) by adopting for the Dark Matter 
the observationally supported Burkert profile, we compute the total angular momentum of the disk and its 
correlation with the rotation velocity; (iii) we find that the distribution function of the spin parameter A peaks 
at a value of about 0.03, consistent with a no-major-merger scenario for the late evolution of spiral galaxies. 

Subject headings: galaxies: halos - galaxies: spiral - galaxies: formation - galaxies: kinematics and dynamics 



1. INTRODUCTION 

The mechanism of galaxy formation, as currently under- 
stood, involves the cooling and condensation of baryons in- 
side the gravitational potential well provided by the Dark Mat- 
ter (DM); in spirals, a rotationally supported disk is formed, 
whose structure is governed by angular momentum acquired 
through tidal interactions during the precollapse phase. 

Under the assumption of specific angular momentum con- 
servation, that holds when the baryons and the DM are ini- 
tially well mixed, the dynamics of the dark halo is directly 
related to the disk scale length (see Fall & Efstathiou 1980). 
This tight connection between halo dynamics and disk geom- 
etry is quantified by the spin parameter A (Peebles 1969). 

The general procedure for the computation of the angular 
momentum has been described in detail by Mo, Mao & White 
(1998); it relies on 3 basic assumptions: (i) the mass of the 
galactic disk is a universal fraction of the halo's; (ii) the total 
angular momentum of the disk is also a fixed fraction of the 
halo's; (iii) the disk is thin and centrifugally supported, with 
an exponential surface density profile. The theory is applied 
to a Navarro, Frenk & White (1997; NFW) DM potential. 

In the present work, we propose to determine the angular 
momentum and the spin parameter of disk galaxies by mak- 
ing use of the observed matter distribution in spirals, and of 
observed scaling relations between halo and disk properties. 
For this purpose, we adopt a modified set of assumptions: we 
relax (i), and use instead an empirical relation that links the 
disk mass to that of its DM halo (Shankar et al. 2005); we 
retain (ii) and suppose total specific angular momentum con- 
servation during the disk formation, i.e., Jo /Mo = Jh/Mh in 
terms of the disk and halo masses Mr>, Mh and of the related 
total angular momenta Jo, Jh', as to (iii), we still assume that 
the disk is centrifugally supported, stable, and distributed ac- 
cording to an exponential surface density profile, but we also 
take into account the gaseous (HI+He) component. Finally, 
we perform the computation for a Burkert halo. 

We adopt a flat cosmology with matter density parameter 
VL M w 0.27 and Hubble constant h = H /\00 km s" 1 Mpc -1 = 
0.71. Given a halo with mass Mh we determine its ra- 
dius as R H = [3M H Q Z M / 47rp e O M (l+z) 3 A H ] 1/3 ; here p e = 
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2.8 x 10 u /i-M Q Mpc is the critical density, Q : M = fW(l + 
z) 3 /[(l-n M ) + n M (l+z) 3 ] and A H = 18 tt 2 + 82(fi^ - 1) - 
39(il z M - l) 2 are the density parameter and the density con- 
trast at redshift z; A H ~ 100 holds at z = 0. 

2. THE ANGULAR MOMENTUM 

The fundamental parameters of the stellar disk and halo 
mass distributions can be obtained straightforwardly by 
means of three observational scaling relations, linking the 
disk mass to the halo mass, to the halo central density, and 
to the disk scale length. 

The total mass of the stellar disk Mr, that resides in a halo 
of mass Mh is given by the relation (Shankar et al. 2005): 
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this holds for halo masses between 10 11 and about 3 x 
10 12 Mq, wide enough to include most of the spiral popula- 
tion, except dwarfs. This relation has been derived by the sta- 
tistical comparison of the galactic halo mass function with the 
stellar mass function; the related uncertainty is around 20%, 
mainly due to the mass to light ratio used to derive the stellar 
mass function from the galaxy luminosity function. 

We model the stellar disk with a thin, exponential surface 
density profile of the form 
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The characteristic scale-length Rd is estimated through: 
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this relation is inferred from dynamical mass determinations 
by Persic et al. (1996). These scale lengths are consistent 
with the data by Dale et al. (1999), Simard et al. (1999), and 
Courteau et al. (2003). 

For the DM, we adopt a Burkert distribution pn(f) = 
Po ^o/( r + ^o)('" 2 + ^o)' w i m the core radius and p {) the 
effective core density. Correspondingly, the total halo mass 
inside the radius r is given by 
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FIG. 1 . — Left panel: the specific angular momentum of the disk as a function of the rotation velocity at 2.2 Rq. Solid line is the result from this work, adopting 
the Burkert profile; dashed line is the best-fit relation from the data collected by Navarro & Steinmetz (2000), see their Figure 3. Right panel: the Tully-Fisher 
relation. Solid line represents the result from this work and dashed line illustrates the fit to the data by Giovanelli et al. (1997). 



with Mo = 1.6/0o/?o being the mass contained inside the ra- 
dius Rq. The density po is determined from the disk mass 
through the relation obtained from the Universal Rotation 
Curve (Burkert & Salucci 2000): 
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For each given virial mass M# (corresponding to a radius Rh) 
we find the density po through Eqs. Q and (|5}; then we nu- 
merically compute the core radius Rq by requiringthat the 
mass Mh(< Rh) inside Rh given by r.h.s. of Eq. J^J equals 
the virial mass Mh- The resulting relation Ro vs. Mh is ap- 
proximated within a few percents by the relation 

log(fl /kpc) « 0.66 + 0.58 log(M ff /lO n M ) ; (6) 

such values of Rq obtained by assuming a mass model out to 
Rh are consistent with those estimated by Burkert & Salucci 
(2000) from the decomposition of the inner rotation curves. 

The total circular velocity of the disk system is V c 2 (r) = 
Vp(r) + V H (r). For a thin, centrifugally supported disk the 
circular velocity is given by Vg(r) = (GM D /2R D ) x 2 B(x/2); 
here x = r/Ro and the quantity B = IqKq- I\K\ is a com- 
bination of the modified Bessel functions that accounts for 
the disk asphericity. Moreover, the halo circular velocity 
is simply V H (r) = GMh(< r)/r, and it is useful to define 
V H = ^GMh/Rh- Given the relations Q, Q, (0 and © 
linking the basic quantities of the system, the shape and am- 
plitude of the velocity profile depend only on the halo mass. 
Note that the uncertainties on these relations combine to give 
a 10% -20% uncertainty on the determination of the velocity 
profile (see Tonini & Salucci 2004). 

In order to check our mass model and empirical scal- 
ing relations, we compute the I-band Tully-Fisher relation 
at r = 3Rd- We obtain the B-band luminosity from the 
stellar disk mass through the relation log(Ls/LQ) w 1.33 + 
0.83 log(M D /M Q ) by Shankar et al. (2005), then convert the 
related magnitude in I-band through the mean colour B—I w 2 
(Fukugita et al. 1995). In Figure 1 (right) we compare the 
result with the data by Giovanelli et al. (1997), finding an 
excellent agreement. 

We compute the angular momentum of the disk as 



Jn = 2-7T 
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with f R = Jo°°x 2 e x V c (xR D )/V H dx, x = r/R D and M D = 
2-kY>qR d . Note that Jo depends linearly both on the mass and 
on the radial extension of the disk, while the DM distribution 
enters the computation through the integrated velocity profile, 
encased into the shape factor /«; the latter slowly varies (by a 
factor 1 .3 at most) throughout our range of halo masses. 

In Figure 1 (left) we show the specific angular momentum 
of the disk vs. the total circular velocity at r = 2.2Ro, com- 
puted as jo = Jd/Mq from Eq. Q. Plotted for comparison is 
also the best-fit relation by Navarro & Steinmetz (2000) from 
their collection of data; note that these authors adopted a flat 
rotation curve, so that /« = 2 and jo =2RdVh- 

We derive the halo angular momentum by assuming the 
conservation of the total specific angular momentum between 
DM and baryons: 

Jh = Jo^ d , (8) 

an ansatz widely supported/adopted in the literature (Mo et al. 
1998; van den Bosch et al. 2001, 2002; Burkert & D' Onghia 
2004; Peirani et al. 2004). Note that small variations of J D 
are magnified by a factor Mh /Mo in the value of Jh, i.e., the 
latter is rather sensitive to the radial extension of the baryons. 

We now consider, along with the stars, the gaseous compo- 
nent that envelops the disk of spiral galaxies. We derive the 
total mass of the gas component from the disk luminosity (see 
above) through the relation 
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by Persic & Salucci (1999), where we have included a fac- 
tor 1.33 to account for the He abundance. Since the gas 
mass is on average much less than the halo mass, both the 
total mass and the rotation curve remain virtually unaltered 
(Vgas ~ -^Mgas/^gas) by the presence of the gaseous compo- 
nent. However, the gas is much more diffuse than the stars, 
reaching out several disk scale lengths (Corbelli & Salucci 
2000; Dame 1993); since most of the angular momentum 
comes from material at large radial distances (van den Bosch 
et al. 2001), we expect the gas to add a significant contri- 
bution to the total angular momentum (Eq. Q), especially in 
small spirals where the gas to baryon fraction is close to 50%. 
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The detailed density profile of the gas in spirals is still un- 
der debate in the literature. However, we are confident that 
the main factors entering the computation of the gas angular 
momentum / gas are just the gas total mass M gas and the ra- 
dial extension of its distribution; in other words, we expect 
that the details of the gas profile do not significantly affect 
the results. In order to check this statement, we considered 
3 different gas models, i.e. (i) a disk-like distribution (DL), 
with scale length aRo, (ii) a uniform distribution (U) out to a 
radius /3Rd', and (iii) anM33-like gaussian distribution (M33; 
Corbelli & Salucci 2000): 



Bullock et al. (2001) proposed the alternative definition 
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where 9 in the second equation is the Heaviside step function. 
As fiducial values of the parameters, we adopt a w 3 in the 
first expression, (3 w 6 in the second one (Dame 1993), and 
ki w 11 .9, k 2 w 5.87 in the last one (Corbelli & Salucci 2000). 
Each profile has been normalized to the total gas mass M gas as 
computed from Eq. l|9). 
In parallel with Eq. |7}, the gas angular momentum will be 

POO 

J m = 2n ^{f)rV c {f)rdr = M m R D V H f m , (11) 
Jo 

where the shape parameter / eas encodes the specific gas distri- 
bution. On comparing its values for the three models we find 
differences of less than 15%, and so confidently choose the 
gaussian profile as a baseline. 

We then compute the halo angular momentum as a function 
of the total baryonic one: 



Jh - (JD + Jgns) 



M D +M R . 



(12) 



The gas is dynamically affecting the system mainly through 
its different spatial distribution with respect to that of the stars, 
adding an angular momentum component that is significant at 
large radii compared to Rp. 

Note that we do not include a bulge component, since it 
would contribute a negligible angular momentum and a mass 
of 0.2Mb at most; in any case, this would slightly lower/// af- 
ter Eq. ill 2b and, as will be evident in the next Section, would 
lower the spin parameter and strengthen our conclusions. 

3. THE SPIN PARAMETER 

The spin parameter is a powerful tool to investigate galaxy 
formation, as it is strictly related to both the dynamics and the 
geometry of the system. We compute its values and distribu- 
tion function based on the results of Section 2. 

The spin parameter is defined as follows: 



A = 
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where G is the gravitational constant, and Eh is the to- 
tal energy of the halo. The latter is computed as \Eg\ = 
2tt J drr 2 pn{r)V 2 after the virial theorem, and having sup- 
posed that all the DM particles orbit on circular tracks. 



./h 



Jh + Jd + Js. 



V2M H R H V H V2(M H +M D +M gas )R H V H ' 



(14) 



where the second equality holds after Eq. d!2l >. We find that 
for Burkert halos the ratio A/ A' is between 1.1 — 1.3 in the 
mass range 10 11 - 3 x 1O 12 M0. In Figure 2 (top panels) we 
plot both A and A' as a function of the halo mass. We highlight 
the difference in the value of the spin parameter when the gas 
component is included, especially in low mass halos. 

To compute the probability distributions V(X) and V(X') of 
the spin parameters, we make use of the galactic halo mass 
function, i.e., the number density of halos with mass Mh con- 
taining a single baryonic core, as derived by Shankar et al. 
(2005). A good fit is provided by the Schechter function 
^(M H ) = (^> /M)(M H /M) a exp(-M H /M) , with parameters 
a = -1.84, M= 1.12 x 1O 13 M and * = 3.1 x 10" 4 Mpc" 3 ; 
note that within our range of halo masses, this is mostly con- 
tributed by spirals. For the computation of "P(A) or "P(A'), 
we randomly picked up a large sample of masses distributed 
according to ^(M h), then compute A or A' for each using 
Eqs. dl3t and dl4> . and eventually build up the statistical dis- 
tributions. During this procedure we have convolved the re- 
lations Jl 31 and dl4> with a gaussian scatter of 0.15 dex that 
takes into account the statistical uncertainties in the empirical 
scaling laws we adopt; these are mostly due to the determina- 
tion of Rd through Eq. 0, for which we have determined the 
scatter by using the disk mass estimates of individual spirals 
reported in Persic & Salucci (1990). 

As shown in Figure 2 (bottom panels), we find a distribution 
peaked around a value of about 0.03 for A and about 0.025 for 
A', when the gas is considered. We stress that this value of A' 
is close to the result of the simulations by D'Onghia & Burk- 
ert (2004), who on average find A' = 0.023 for spirals quietly 
evolving (i.e., experiencing no major mergers) since z ~ 3, see 
their Figure 4. In addition, Burkert & D'Onghia (2005) argue 
that this value of A' provides a very good fit to the observed 
relation between the disk scale length and the maximum rota- 
tion velocity (see their Figure 1). We also stress that our most 
probable value for A is in agreement with the results by Gard- 
ner (2001), Vitvitska et al. (2002) and Peirani et al. (2004), 
who find a peak value of the distribution function at around 
0.03 for halos that evolved mainly through smooth accretion. 

4. DISCUSSION AND CONCLUSIONS 

In this Letter we have computed the angular momentum, 
the spin parameter and the related distribution function for 
DM halos hosting a spiral galaxy. We have relied on observed 
scaling relations linking the properties of the baryons to those 
of their host halos, and have assumed the same total specific 
angular momentum for the DM and the baryons. 

Our main findings are: (i) we show that including the gas 
component beside the stars has a remarkable impact on the 
total angular momentum; (ii) by adopting for the DM the ob- 
servationally supported Burkert profile, we compute the total 
angular momentum of the disk and its relationship with the 
rotation velocity; (iii) we obtain A' m 0.025 and A w 0.03 as 
most-probable values of the spin parameters. 

Simulations based on the ACDM framework, performed by 
various authors (Bullock et al. 2001; D'Onghia & Burkert 
2004), have shown that the distribution of the spin parameter 
A' for the whole halo catalogue peaks at around 0.035, signif- 
icantly higher than our empirical value. However, D'Onghia 
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FIG. 2. — The spin parameter and its distribution function. Top panels: A' (left) and A (right) as a function of the halo mass, when the gas component is included 
in the system (solid line) and when it is not (dashed line). Bottom panels: the distribution function of A' (left) and A (right), again with gas and without gas. 



& Burkert (2004) highlight that if one restricts one's atten- 
tion to halos that hosts spirals and have not experienced major 
mergers during the late stages of their evolution (z < 3), the 
average spin parameter A' turns out to be around 0.023, very 
close to our observational result. Moreover, Gardner (2001) 
and Peirani et al. (2004) showed that the spin parameter A 
undergoes different evolutions in halos that have grown up 
mainly through major mergers or smooth accretion: in the for- 
mer case A takes on values around 0.044, while in the latter 
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case A has lower values around 0.03. 

Thus our findings point towards a scenario in which the late 
evolution of spiral galaxies may be characterized by a rela- 
tively poor history of major merging events. 
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